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Abstract 

Recently, a finite-temperature real-time static potential has been introduced via a 
Schrodinger-type equation satisfied by a certain heavy quarkonium Green's function. Fur- 
thermore, it has been pointed out that it possesses an imaginary part, which induces a finite 
width for the tip of the quarkonium peak in the thermal dilepton production rate. The imagi- 
nary part originates from Landau-damping of low-frequency gauge fields, which are essentially 
classical due to their high occupation number. Here we show how the imaginary part can be 
measured with classical lattice gauge theory simulations, accounting non-perturbatively for 
the infrared sector of finite-temperature field theory. We demonstrate that a non-vanishing 
imaginary part indeed exists non-perturbatively; and that its value agrees semi-quantitatively 
with that predicted by Hard Loop resummed perturbation theory. 
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1. Introduction 



The notion of a static potential, generalizing the potential that appears in the Schrodinger- 
equation of non-relativistic quantum mechanics, is thought to play a role for heavy quarko- 
nium physics in QCD. If the energy of a two-quark system, E, is close to twice the mass of the 
heavy quark, M, so that the combined "kinetic energy" of the two quarks, E — 2M, is small 
compared with M, then we may assume the quarks to be "static" to a good approximation, 
moving only slowly in the attractive potential generated by the colour fields. 

To turn this intuitive picture into a quantitative description requires the use of effective field 
theory methods. At zero temperature various energy and momentum scales can be identified, 
the small expansion parameter being related to the ratio {E — 2M)/M. The relevant effective 
theory is called NRQCD [Tj, or one of its variants, like pNRQCD [2]; for reviews on the various 
effective theories used for describing heavy quarkonium, see refs. [3ll3|. The static potential 
plays the role of a certain matching coefficient in these effective theories: it is related to, but 
not identical with, the non-perturbative static potential that is traditionally defined from a 
large Euclidean Wilson loop in lattice QCD. 

At finite temperatures, the situation becomes more complicated than at zero temperature. 
Indeed, finite-temperature field theory possesses many momentum and energy scales of its 
own: gluonic momenta could parametrically be /c ~ irT, gT, g'^T [5l[6], while gluonic frequen- 
cies (energies) can be even softer, down to ~ g^T [TJ [8]. Here T is the temperature and g 
is the QCD gauge coupling. The relevant effective description now depends on the relation 
of these scales to the scales already appearing in the zero-temperature situation. 

In fact, at finite temperatures, the situation is quite complicated even at the leading non- 
trivial order in g. This might be anticipated from the fact alone that the definition of a static 
potential based on the Euclidean Wilson loop appears to lose its meaning: the Euclidean time 
direction becomes compact, and large Wilson loops do not possess the same interpretation 
as at zero temperature. Replacing the Wilson loop by a correlator of Polyakov loops does 
not remedy the situation [9l|10]. Moreover, physics lives in Minkowski spacetime, which at 
finite temperatures in general requires a non-trivial analytic continuation 

Recently an attempt was made to give a proper definition of a static potential in this situa- 
tion, in the sense of obtaining an object which has a direct connection to the spectral function 
of the heavy quarkonium system (at least up to some order in perturbation theory) [12]. For- 
mally, the static potential could be defined as a certain coefficient in the large-M expansion 
of an equation of motion satisfied by a suitable heavy quarkonium Green's function. At 
leading non-trivial order, the corresponding object was computed in Hard Thermal Loop [13] 
resummed perturbation theory in ref. [12]. It was found that, at least to this order, the static 
potential can also be obtained from a specific analytic continuation of the Wilson loop de- 
fined in Euclidean spacetime with a compact time direction. At the same time, this analytic 
continuation yields properties that are not familiar from the zero-temperature context: in 
particular, the potential develops an imaginary part. 
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It is the purpose of the present paper to elaborate on the existence of an imaginary part. 
We start, in Sec. [21 by reviewing the definition(s) introduced in ref. [12j. In Sec. [3] we 
argue that the imaginary part of the static potential remains non-zero in the classical limit, 
by computing it perturbatively in classical lattice gauge theory. Given that perturbative 
computations at finite temperatures may ultimately suffer from infrared divergences, we carry 
out non-perturbative Monte Carlo simulations in classical lattice gauge theory in Sec. HI and 
compare the results with those of the perturbative computation. We conclude in Sec. El 



2. Definition of a real-time static potential 

We start by defining a certain Green's function in hot QCD. Let r be a point-splitting vector, 
and a generic heavy quark field operator in the Heisenberg picture. Then we introduce 

(7>(t,r)^ y"d3x(|(t,x + 07'^Ty^(t,x-0 |(0,0)7;.V^(0,0)) , (2.1) 

where is a Wilson line along a straight path connecting the adjacent operators, inserted in 

order to keep the Green's function gauge- invariant; the metric is r]fj,iy = diag(H ); and the 

expectation value refers to (...) = Z~^Tr: [ex.p{— H /T) {...)], where Z is the partition function, 
H is the QCD Hamiltonian operator, and T is the temperature. 

The significance of the Green's function in Eq. (|2.1|) is that if we take the limit r ^ 0, and 
subsequently Fourier transform with respect to the time t, then we obtain a function which 
is trivially related to the heavy quarkonium spectral function, p{uj), in the vector channel: 

p(a;) = i(l-e-^) y""dte-*C'>(t,0). (2.2) 

On the other hand, keeping r for the moment, makes it easier to analyse this Green's 
function in perturbation theory. 

Let us consider C'>(t, r) in the limit that the heavy quark mass M is very large. Then 
C'>(t,r) satisfies a Schrodinger equation of the type 



idt - 

with the initial condition 



V2 1 



2M + VAt,r)-j^ + Oi^ 



C'>(i,r) = 0, (2.3) 



(7>(0,r) = -6iVcJ(3)(r) + o(l:) . (2.4) 



The terms shown explicitly in Eqs. (j2.3p . ()2.4p result from a tree- level computation; they also 
develop multiplicative radiative corrections which we have omitted for simplicity. In contrast, 
the potential denoted by V>(t, r) originates only at 1-loop order. It can be defined as the 
coefficient scaling as 0{M^), after acting on C'>(t, r) with the time derivative idt- 
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Now, as Eq. (j2.3p shows, V>(t, r) can even be defined in the limit M — > oo, provided that 
the trivial factor 2M is shifted away by a redefinition of time, as is standard in NRQCD. 
In ref. [12j, the computation in this limit was carried out to 1-loop order, 0{g'^), in Hard 
Thermal Loop resummed perturbation theory [13j. It was found that at this order Vy{t,r) 
can in fact be extracted from the equation 



idtCE{it,r) = Vy{t,r)CE{it,r) 



(2.5) 



where the function CE{T,r) is nothing but the Euclidean Wilson loop, computed with an 
imaginary time coordinate r, with gauge fields periodic in r ^ r + h/T. 

The expression that was obtained for Vy{t,r) in ref. [T2j reads (the superscript refers to 
the order in g; we keep 1; and we assume the use of dimensional regularization) 



dp° 

TT 



mo + 
d3p 



exp(— TTinr) 



(2.6) 



(2vr)^ 



(1 — cospsr) X 



~i|p"|t 



+ nB(b°|) e 



-i\p°\t 



(2.7) 



Here Cp = (-/V^ — l)/2A'c; mj^ is the Debye mass parameter (actually of dimensionality 
1/distance rather than mass); n-Q[x) = l/[exp(^x/r) — 1] is the Bose distribution function; 
and we have chosen r = (0, 0,r). The r-independent term in Eq. ()2.6p amounts to twice 
a thermal mass correction for the heavy quark. For the gluon spectral functions pe^Pt we 
assume the conventions specified in appendix A of ref. |14j. 

It can be observed that Eqs. (j2.6p . ()2.7p contain both a real and an imaginary part. In 
particular, the familiar- looking structure in Eq. (|2.6p . representing a Debye-screened Coulomb 
potential, is real, while the manifestly "thermal" part in Eq. (j2.7p . containing ne, is purely 
imaginary. As pointed out in ref. jl2], this purely imaginary term remains non-zero in the 
limit t — > oo, because of Bose-enhancement at small frequencies, h\p^\ <^ T (cf. Eq. (12. 9p ). 

The Bose-enhanced term corresponds to the physics of the classical limit /i — > 0, in which 
situation riedp'^l) = T/h\p^\. In fact, all other terms vanish in this limit, being multiplied by 
h. The classical potential then reads 



g^CpT 



d^p 
(2^ 



(1 — cospsrj / 

vr 



lim 



1 



1 



pe(/,p) + 



1 



1 



Pt(p°,p) 



p2 (p0)2^ 

where we have also simplified the way in which p'^'s appear. For large times. 



lim 

i— >oo 



P' 



,0 



2TTi (5(p°) 



(2. 



(2.9) 
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and we obtain 



Ff(^,r) = 2z5^C^r|^(l-cosp3r)^limJim^^^^ (2.10) 

If we also take the limit r — > oo, the cosine-term in Eq. (j2.10p drops out. Assuming for a 
moment that the two limits and the integration in Eq. (j2.10p commute, and making use of 
the known Hard Thermal Loop form of Pe{p^tP) (Eq. (B.13) of ref. [12] shows pE at small 
\p^\ with our conventions), containing the parameter m^, then leads to the provisional result 



J (27r)'^ IpKp^ + m^)^ h^o Att Ait 



Note that in the quantum theory the same result is obtained for the asymptotic value 
6V^ '{oo, oo) [H], and in this sense Eq. (j2.1ip is indeed the correct physical expression. 

Now, given that the integral in Eq. (j2.1ip is finite, it might be assumed that the result is 
independent of the regularization procedure. It turns out that this argument is too naive: in 
fact, diverges as g^T"^ /h in dimensional regularization, indicating that the classical limit 
may introduce ultraviolet singularities. In particular, if the ultraviolet is regularized by a lat- 
tice rather than dimensionally, with a spatial lattice spacing a, then the limit becomes finite, 
lim^i^o "^-D °^ g^T/a [EJ [E]. Thus the classical limit in Eq. ()2.10p does exist, but the price 

(2') 

to pay is that pE and subsequently V^^ (oo,cx3) depend on the details of the regularization 
procedure. In particular, carrying out the limits in the order indicated by Eq. (|2.10p with 
lattice regularization, does not lead to the expression in Eq. ()2.1ip (cf. Fig. [3] below) . 

Fortunately, this problem is not too serious: an analogous situation was met in studies of 
the sphaleron rate in the electroweak theory, yet classical lattice gauge theory simulations [T7] 
did yield non-perturbative physical information, once properly interpreted (see, e.g., refs. [181 
[19]). In our case. Fig. [3] implies that we cannot use classical lattice gauge theory simulations 
to compute corrections directly to Eq. (12. lip . However, we can compute the analogues of 
Eqs. ()2.8p - (j2.1ip with Hard Thermal Loop perturbation theory adapted to the ultraviolet 
physics of the classical lattice [ISj [l6] , and compare subsequently these results (the dashed 
curve in Fig. [3]) with a non-perturbative determination. In this way we can probe the infrared 
sector of thermal field theory, which indeed is classical in nature. 



3. Perturbative real-time static potential in classical lattice gauge theory 

We assume that the theory is regularized by introducing a cubic spatial lattice, while the time 
coordinate is continuous. Gauge field configurations are generated with a Wilson-discretised 
Hamiltonian (cf. Eq. ()4.ip below), and evolved with the classical equations of motion (cf. 
Eqs. (I4.3p . ()4.4p below). The results depend on a single parameter, 

(3.1) 
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Figure 1: Integration contours for the classical real-time static potential. 



where Ca = N^. 

The way to carry out perturbation theory in this situation was worked out in refs. [151 116] . 
The procedure is analogous to Hard Thermal Loop resummed perturbation theory [13], with 
technical differences originating from the different ultraviolet physics. We will refer to this 
procedure as Hard Classical Loop (HCL) perturbation theory. 

Let us start by introducing the notation 



Pi = - sm I 

a 



i=l 



i=l 



Also, the integration measure is denoted by 

dp = 



Pi 



(3.2) 



ir/a ^3 



-7r/a(2vr)^ 



(3.3) 



Then, we consider Eq. (12. Sp . with a few straightforward modifications following from the 
introduction of lattice regularization: 













(p°)0 



dp" 



Here, 



Pe{p^ p) = wr ^e{p^ + ie, p) - Ae{p^ - ie, p) 
2^ L 



with e = 0"*", and the propagator A^; has the form 

Ai^(/,p) 



1 



p2_(pO)2 + nE(pO,p) ■ 

The limit lim/^^o is assumed everywhere but not shown explicitly. 



(3.4) 
(3.5) 

(3.6) 



3.1. Behaviour at finite times 

In order to simplify Eq. (j3.4p . it is convenient, following ref. [20], to view the p'^-integration 
as an integral in the complex plane, and to deform the contour suitably. As it stands, 
the integrand in Eq. (|3.4p is finite at p*^ = (note that p^; is linear in around the origin). 
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However, the part multiplying pE contains a pole at = 0; this pole just does not contribute 
because of the mentioned property oi pE- Writing p^; as in Eq. (j3.5p . this means that we can 
view the original integral as indicated in the left-most drawing in Fig. [H and then also deform 
it accordingly. Subsequently, e can be taken to be finite, because there are no singularities 
outside of the real axis. Furthermore, assuming t > 0, terms multiplied by exp + ie)t] 
and exp[— — ie)t] must vanish, because we can imagine taking e arbitrarily large. The 
integrand can only decrease in this limit, and being multiplied by exp(— ei), the integral then 
vanishes (it is a good cross-check of the numerics to verify the vanishing at any finite e). 
Finally, the symmetry properties of the integrand allow to reflect the lower of the remaining 
contours to the upper half-plane. We thus obtain 



Im 



2Cf 



dp(l 
1 



cosp3r)<^ tAE{0,p) 



d/ i(pO+i,)i 



As(/ + ie,p) + 



Pi 



— oo 

1 



27r 



AT(p° + ie,p) 



(3.7) 



where the first term is the contribution of the pole in Fig. [TJ 

Let us stress that the integration in Eq. (|3.7p is independent of the value of e > 0, since 
there are no poles in the upper half-plane. Checking the independence in practice offers 
another cross-check for the accuracy of the numerical integration. Naturally, small values of 
e are difficult, because the integrand becomes strongly peaked around the origin, while large 
values of e are also difficult, because the latter term is multiplied by exp(et), whereby the 
numerical errors of the integration are exponentially amplified at large t; a useful compromise 
appears to be e ~ 1/a. In general, it is advantageous to decrease e when increasing t. 

In order to insert the propagators A^;, At, we need to know the self-energies II^;, IIj' (cf. 
Eq. ()3.6p ). Starting from the spatial part of the gluon self-energy [151 [16] 



Uij{p'>,-p) = 2g^TCA /dq 



1 



q2 pO 



p V 



with 



(3., 



(3.9) 



and employing the projection operators Pj^u^Pj^u defining He^^t (we use the conventions 
specified in appendix B of ref. fl2]), we obtain 

(/)2\ / S 



ni^(p°,p) 
nT(/,p) 

Here 



2g'TCA[ 1 



P" 



A'Ka 



p 



,0 



(^T 

p2 



Ana 



p 



,0 



Ana 



q^ — p • V 



dq-4 



p • V, 
+ / dq 



q 



p- V 



(3.10) 
(3.11) 

(3.12) 
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where S ~ 3.175911535625 is a trigonometric factor which can be expressed in terms of the 
complete elhptic integral of the first kind [21] . Note that with the form in Eq. (j3.10p . the 
combination containing A^; in Eq. (13.71) becomes 



1 



p2 (pO + ie)2 



As(p°+ie,p) 



1 



(p° + ie) 



(3.13) 



Denoting the square brackets in Eq. (j3.7p by T{p^ + ie, p), and making use of the properties 



T(-pO + ie, p) = J(/ - ie, p) = T(p° + ie, p) 



(3.14) 



the integral over can furthermore be reduced to regular cosine and sine transforms: 



27r 

et /-oo 



e ' "-'"l{jP + ie, p) 



ei /-oo 

— y dp° |cos(p°t) Re p(p° + ie, p)J + sin(p°t) Imp(p° + ie, p)J | . (3.15) 



Though efficient routines for such transforms exist, it is also clear that the accuracy require- 
ments grow exponentially with t, so that very large times are difficult to reach. 

In a practical lattice study, the system possesses not only a finite lattice spacing, but also 
a finite extent, L = Na, where is the number of lattice points. We assume that the box is 
cubic and that periodic boundary conditions are imposed in every direction. Furthermore, let 
us assume that we use changes of integration variables to write the momentum integrations 
over the "positive" octant only, 



IX (a 



27r 



-^Q{Pi), g{p^)^ HPi)+H-p^) 



1,2,3 



(3.16) 



In a finite volume this then goes over into 



"/'^dK 1 



27r 



Na 



N/2-1 

1 = 1 



a' 



(3.17) 



Note that in finite volume, the analytically known integral in Eq. ()3.12p should also be 
replaced by a numerically evaluated sum. 

Now, Eq. (I3.17|) contains also a contribution from the zero-mode, pi = 0. Its treatment 
requires in general some care. In 11^, 11^, loop momenta are by definition "hard" : it is sensible 
(and in fact necessary) to leave out the zero-mode. In the remaining sum in Eq. (j3.7p . in 
contrast, momenta could be soft: we thus keep the contribution of the zero-mode as well 
(even though the practical effect is small). Note that for the zero- mode. 



n£;(/,o) = nT(p°,o) 



ijj. 



pi 



pi 



-g'TCA /dq 



q 



2^2 



(3.18) 
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|3= 16, N= 12, analytic 

' \ ' \ ' r 




r/a= 1 

-- r/a = 2 

■■■■ r/a = 3 

— r/a = 4 



3 6 9 12 

t/a 



Figure 2: The imaginary part of the classical real-time static potential, to leading non-trivial order 
in HCL-resummed perturbation theory (Eq. (|3.7p with /dp replaced by a finite- volume sum), for 
(3=16, N = 12, Nc = 3. 

where the integration can be replaced by a sum (without zero-mode) as before0 Then the 
combination in Eq. (j3.7p becomes 

(1 — cos p^r) 



1 1 \ 1 /II 



p2 (pO ie)2 y p2 _ (pO _^ ^g)2 _^ ^2^ \^^2 p2 ^ p2 _ (^0 _^ ^g-)2 _^ ^2^ 



p=o 1 



' 2 -(p0 + ie)2+^2^ • (3.19) 

In practice, for the values > 12 that we have used, finite-volume effects are almost 
invisible at small times. Perturbative finite-volume effects grow rapidly with time, however, 
and also with distance. On the other hand, perturbation theory tends to overestimate their 
significance, since it lacks the mass gap generated by the confining dynamics. At the same 
time, whenever justified, it appears to be numerically advantageous to use the finite-volume 
expressions, which contain a six-fold exact sum, rather than to approximate the corresponding 
infinite-volume continuous six-dimensional momentum integration numerically. Therefore we 
plot the perturbative expression only in the range where the perturbative finite- volume effects 
are small, t/a< 10; an example of a result is shown in Fig. [2l 



^In infinite volume, tJpj = g^TCA(3E/27r - l)/6a, where E is the constant in Eq. ((3TT2|) [20] 
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3.2. Value in the large-time limit 



As mentioned, it is not easy to evaluate numerically the HCL-resummed perturbative expres- 
sion in Eq. ()3.7p . once the time coordinate becomes large: e should be decreased, whereby 
the integrand becomes strongly peaked; and one should replace the finite-volume sums with 
infinite- volume momentum integrals, whereby the numerical cost increases. To get a handle 
on this limit we can, however, proceed in another way, without making use of the contour 
trick, and thereby obtain the correct version of Eq. (j2.1ip on an infinite spatial lattice. The 
starting point is then Eq. (j2.10p . 

For infinitesimally small e, the small-p'^ behaviour of II^; reads (cf. Eq. ()3.10p ) 



+ inp 



0\2\ 



(3.20) 



Making use of the definition in Eq. ()3.5p . the asymptotic value from Eq. ()2.10p can be written 
as 



Im 



(2), 



cl 



oo, r 



ttCfC\ 



d^x 



1 — cos(7rx3r/a) 



1 ^ (y2)V2 



(3.21) 



where we have gone over to a notation where the integration variables are made dimensionless 
by going to lattice units, and the integration range is restricted to the unit box and its 
reflections: 



2 sin 



Xi 



sin(7rxi) 



Xi G 



-1,1). 



(3.22) 



Moreover, we have made use of the symmetry of the integrand, in order to restrict the 
integration to positive Xj. 

Among the eight octants of the y-integration, the (5-function gets realised in six only, and 
by changes of integration variables we can combine all the contributions together: 



Im 



Kf^(oo,0 



2ttCfCI 



f3 



d'^x 



3 — cos(7rxir/a) — cos(7r2;2r/a) — cos(7rx3r/a) 







(i2 + C72SA/3) 



3 S{xiyi+X2y2 
d y — 



xsy-i) 



(y2)l/2 



(3.23) 



It is now straightforward to carry out the integration over, say, X3, to remove the (5-function, 
and also to make use of the symmetry of the remaining integrand in xi ^ X2- The 5- 
dimensional integral left over can be evaluated numerically without too much trouble. The 
result is shown in Fig. [3l for A'c = 3 and the case r/a = 00, when the cosine-term does 
not contribute. Values of Eq. (j3.23p at finite r/a can be found in Table 1 below. Comparing 
Table 1 with Fig. [3l we observe that for, say, /? = 16, the distance r/a = 4 gives a value which 

is already relatively close to the asymptotic one. To summarise, in resummed perturbation 

(2) 

theory Im[l/[ (00, r)] is definitely non-zero at all r 7^ 0. 
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Figure 3: The asymptotic value of the real-time static potential Vc\{t,r), on an infinite lattice, 
to leading non-trivial order in HCL-resummed perturbation theory (Eq. (I3.23p ). for = 3. For 
comparison, we also show the expression on the right-hand side of Eq. (|2.1ip , which corresponds to 
the asymptotic value of the potential {t, r) in the full continuum quantum theory. 

4. Classical lattice gauge theory simulations 

The computation in the previous section was based on resummed perturbation theory, but 
it was only carried out to a fixed order. Let us try to estimate the expansion parameters of 
such a computation. Using continuum notation, the vertices of each new loop order bring in 
a factor g^T. At the same time, the mass scales appearing in the dynamics are the ultraviolet 
cutoff scale, A 1 /a, as well as the confinement scale of three-dimensional Yang-Mills theory, 
^ g^T [3 El- The resummation accounts for the dominant influence of the hard cutoff 
scale A on the dynamics of the infrared modes, and is associated with corrections of the 
type g^T/K ~ 1//?. However, it does not account for the self- interactions of the infrared 
modes, which may lead to an expansion parameter of the type g'^T/niQ, ~ 1. Therefore, we 
would like to compare the resummed perturbative result with a non-perturbative numerical 
computation. 

The practical procedure of generating classical gauge field configurations is the follow- 
ing [T7|. Since our observable will be gauge-invariant, we may choose a gauge; it is convenient 
to work in the temporal gauge, C/j(x, t) = 1. Here C/^ is an SU(A'c) link matrix. The canoni- 
cal degrees of freedom are the spatial link variables f7j(x, t) and the matrix- valued "canonical 
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momenta" Ei{x,t), which transform in the adjoint representation. We denote the generators 
of the gauge group by T'*, and assume them normahsed through Tr [T°'T^] = 6"'^/2. Further- 
more, Ei = EfT'', Sij{x) = Uj{x)Ui{x + j)uj{x + i), where i = aci, and C/_i(x) = C//(x-i). 

With this notation, the procedure starts by generating initial configurations (at time t = 0) 
according to the partition function |l7j 

Z = JvU,VE,6{G)expL(3Y, Y.{^ - -^^^^^ Pij) +Y.^HEf) \, (4.1) 

X '-i<j i 

where Pij is the spatial plaquette, and the Gauss law function reads 

G(x,t) = J^[s,(x,t)-[/_i(x,t)S,(x-i,t)[/l.(x,t)] . (4.2) 

i 

To obtain configurations extending to t > 0, we solve the equations of motion 

adtUiix,t) = i{2CA)^Ei{x,t)Ui{x,t) , (4.3) 

1 

adtEl{^,t) = -(-^y ImTT[T'U,{^,t)J2sU^,t)] . (4.4) 

These four-dimensional configurations are then used for evaluating the real-time observables. 
In all that follows, we fix Ca = Nc = 3, even though we have also carried out some simulations 
at iVc = 2 crosscheck. 

It is worth stressing that in Eqs. (j4.3p . (j4.4p . the lattice spacing is finite in spatial directions 
only. In practice, of course, the time direction needs to be discretised as well, but with a 
very small lattice spacing, at <C a. As a check of the time evolution, it is useful to control 
the conservation of the Gauss law and of the total energy. 

To specify the observable to measure, we adopt the definition in Eq. (|2.5p as our non- 
perturbative starting point. The object appearing here is a specific analytic continuation of 
the Euclidean Wilson loop, and corresponds formally to a time ordering generally denoted 
with the subscript (...)> [22] : 

Cy{t,r) = CE{it,r) . (4.5) 

At the same time, the classical {h 0) part of the analytic continuation, which we denote 
by Cci{t, r), is independent of time ordering (lim^i^o C> = lima-^o C< = lim^i^o [C'> + C<]/2)- 
In fact, Cci(t, r) is nothing but the classical Wilson loop, defined in Minkowski time. (Note 
that having chosen the gauge Ut = 1, the classical Wilson loop amounts really to a two-point 
correlation function of two spatial Wilson lines, both of which are local in time.) The classical 
static potential is then measured from 

idtC,i{t,r) = V,i{t,r)C,i{t,r) . (4.6) 

It turns out that Cci{t,r) is real for all times (within statistical errors), and slowly decaying. 
Therefore, Vc\{t,r) is purely imaginary, with a negative imaginary part. 



11 



P = 16, N = 12, simulation 




t/a 

Figure 4: The classical Wilson loop, measured with classical lattice gauge theory simulations, as a 
function of time in units of the spatial lattice spacing, for /3 = 16, N — 12, — 3. 

The technical implementation of our simulation follows earlier work [n\ I18j . However, 
to speed up thermalization, we have implemented the idea mentioned in ref. [23], whereby 
the link variables Ui are first pre-thermalized with regular Monte Carlo techniques in the 
dimensionally reduced SU(3) + adjoint Higgs theory (we use the code described in ref. ^^)- 
Since it is non-trivial to match the parameters of that theory and our effective theory exactly, 
those configurations are not yet fully thermalized. However, this is not a problem, they now 
need only to be evolved for a short time a la refs. [IT^ ITSj . in order to reach the correctly 
thermalized configurations corresponding to the exact parameters of Eq. (j4.ip . 

We have carried out simulations mostly with f3 = 16; since analytic HCL predictions also 
refer to a finite value of /?, there is no need to carry out a continuum extrapolation (cf. 
Fig. [3]). As typical lattice extents we have used = 12 and N = 16; the difference of the 
results between these two is only at the percent level (cf. Table 1 below). The time variable 
is discretised with a spacing at, with a value at/a = 0.01; measurements are recorded every 
10th time step. We stress that thermalization is only carried out in the beginning, while the 
subsequent time evolution is deterministic and follows Eqs. (14. 3p . (14. 4p . 

A representative result for the classical Wilson loop is shown in Fig. HI The corresponding 
potential, extracted from Eq. (j4.6p . is shown in Fig. [5l The result can be compared with 
Fig. [21 showing the HCL-resummed perturbative prediction with the same parameter values. 
The general shapes are seen to match each other to a remarkable degree. On closer inspec- 
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P = 16, N = 12, simulation 




3 6 9 12 

t/a 



Figure 5: The imaginary part of the real-time static potential, measured with classical lattice gauge 
theory simulations, as a function of time in spatial lattice units (with the same parameter values as 
in Fig. 3]). Vertical lines indicate statistical errors, but they are almost invisible in this time range. 

tion, however, the amplitude of the oscillations is larger in the simulation; the frequency of 
oscillations is smaller (the oscillation period is larger); and the absolute value of the potential 
is larger (the imaginary part is more negative). 

In order to quantify the difference, we note that at large times, the potential obtains a 
constant value (or, in terms of the Wilson loop, Cd decays exponentially, cf. Fig. We 
estimate this value by fitting a constant to data in the range t/a = 15... 30, where initial 
transients have died out, yet the statistical errors are still relatively small for all parameter 
values that we have used. The results of the fits are shown in Table 1. Unfortunately, 
statistical errors rapidly increase with r/a, and we are not able to go to large enough values 
for the r-dependence to have flattened off. Nevertheless, the values at r/a = 4 already indicate 
that the asymptotic value is larger (in absolute magnitude) than the analytic HCL estimate, 
by some 50 - 100%. The fact that there thus appears to be somewhat more "damping" in 
the non-perturbative classical dynamics than in the HCL estimate is not surprising: other 
observables have yielded indications of a similar pattern [20] . 

We have also carried out so-called HTL simulations, both with the implementation based on 
treating the velocities of the hard particle degrees of freedom with spherical harmonics [19\ |25] , 
and through a discretization based on platonic solids [26]. The HTL-simulations introduce 
a new parameter, which we refer to as m-n"'^'''. For small mn"''\ say ^ 0.2, the 
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parameters 


value of Im[yci(oo, r)/ g^T\ 


/? 


N 


(bare) 


method 


confs 


r/a = 1 


r/a = 2 


r/a = 3 


r/a = 4 


16.0 


12 


0.00 


simulation 


200 


-0.060(2) 


-0.156(8) 


-0.246(26) 


-0.319(56) 


16.0 


16 


0.00 


simulation 


160 


-0.059(2) 


-0.155(8) 


-0.245(22) 


-0.326(48) 


16.0 


12 


0.21 


simulation 


200 


-0.059(2) 


-0.147(7) 


-0.229(23) 


-0.297(51) 


16.0 


12 


0.35 


simulation 


182 


-0.030(2) 


-0.064(5) 


-0.096(12) 


-0.118(21) 


13.5 


12 


0.25 


simulation 


142 


-0.071(2) 


-0.174(10) 


-0.270(33) 


-0.341(97) 


16.0 


oo 


0.00 


analytic 




-0.0601 


-0.1145 


-0.1507 


-0.1737 



Table 1: The asymptotic values Ini[V^i(oo,r)/g^r], obtained by fitting a constant to data in the 
range t/a = 15. ..30. The numbers in parentheses indicate the uncertainties of the last digits. The 
cases aTOQ^^"^"'' — 0.00 refer to the classical theory without HTL degrees of freedom. The bottom row 
gives the perturbative values from Eq. (|3.23p : the perturbative result at r/a = cx) is -0.2152. 

results are practically identical with those of the classical simulations (cf. Table 1). With 
increasing rn-n'"^"', say arn-n'"^"' > 0.35, we see some discrepancies; in particular, the asymptotic 
value Im[yci(oo, r)] decreases in absolute magnitude (cf. Table 1), as one would expect in a 
situation where a lattice-induced dynamical Debye screening is overtaken by a continuum- 
like parameter (cf. Fig. [3l). At the same time, this method is not really bringing us closer to 
the physical continuum limit (i.e. the continuum limit of the quantum theory), since in that 
situation the bare parameter [m^'"^''']^ should in fact become negative as /5 is increased, in 
order to cancel ultraviolet divergences from the dynamics. Unfortunately, the implementation 
of HTL-simulations that we have followed, based on refs. [19\ I26j. does not allow to simulate 
at [mD""^"']^ < 0, and a single bare parameter would in any case not allow to renormalise all 
the observables that can be measured with classical lattice gauge theory |15[ [TB] . Therefore, 
we omit a more detailed discussion of the HTL-simulations from here. 

5. Conclusions 

The purpose of this paper has been to elaborate on the fact that the finite-temperature real- 
time static potential extracted from an analytic continuation of the Euclidean Wilson loop, 
which can (at least to some order in perturbation theory) be inserted into a Schrodinger-type 
equation governing the behaviour of a certain heavy quarkonium Green's function, contains 
an imaginary part. As discussed elsewhere [U], this imaginary part has an impact on the 
heavy quarkonium spectral function at temperatures above a few hundred MeV, introducing 
a width to (the tip of) the resonance peak. 

Physically, the imaginary part implies that quarkonium at high temperatures should not be 
thought of as a stationary state. Rather, the norm of its wave function decays exponentially 
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with time. 

The imaginary part emerges from Bose-enhanced infrared dynamics and, in field-theory 
language, is classical in nature (in particle language, it corresponds to a net disappearance 
of low-energy off-shell gauge particles, due to inelastic 2 — > 1 and 1^2 scatterings with the 
hard particles in the plasma). We have computed the imaginary part with Hard Classical 
Loop resummed perturbation theory, and with non-perturbative classical lattice gauge theory 
simulations. The comparison of the results, Fig. [2] and Fig. 5, or Table 1, shows reasonable 
qualitative agreement. We conclude that non-perturbative colour- magnetic fields do not play 
a dominant role for the imaginary part of the real-time static potential; however, the non- 
perturbative corrections, together with higher-order perturbative terms, are important on 
the quantitative level, bringing about some 50 - 100% increase in the absolute value of the 
imaginary part at large times (for /? = 16), akin in magnitude to the correction observed 
for the static Debye screening mass [24! • case, our study confirms that an imaginary 

part exists, and suggests that Hard Loop perturbation theory presumably provides for a 
reasonable first estimate for it also in the full quantum theory. 

For physical applications, such as determining the quarkonium spectral function, it is essen- 
tial to use the full quantum theory, rather than the classical one. Moreover, it is convenient 
to use dimensional regularization. Finally, as argued in ref. [13], the static potential should be 
evaluated at t ^ r. The perturbative static potential in this limit, (oo, r), including both 
a real and an imaginary part, can be found in Eqs. (4.3), (4.4) of ref. [12j, and has already 
been employed for estimating the quarkonium spectral function in ref. [14]. Increasing the 
imaginary part by some 50 - 100% in the results of ref. [13] lowers and widens the quarkonium 
peak, but the effect is not dramatic; in general, it appears that the spectral function is more 
sensitive to the real than the imaginary part of the real-time static potential. 

As the next step of our program, we would therefore like to get a non-perturbative handle 
also on the real part of the real-time static potential, Vy{t,r), entering the Schrodinger- 
equation. In particular, it would be important to clarify its connection to the other static 
potentials that are being used for studying the spectral function of heavy quarkonium in 
high-temperature QCD (for recent work and references see, e.g., ref. [27j). 
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